Normal and Lateral Casimir Forces between Deformed Plates 



Thorsten Emig 

Institut fiir Theoretische Physik, Universitdt zu Koln, Ziilpicher Strafle 77, D-50937 Koln, Germany 

Andreas Hanke 
Institut fiir Theoretische Physik, Universitdt Stuttgart, 
Pfaffenwaldring 57, D-70550 Stuttgart, Germany 

Ramin Golestanian 

Institute for Advanced Studies in Basic Sciences, Zanjan 45195-159, Iran 
and Institute for Studies in Theoretical Physics and Mathematics, P.O. Box 19395-5531, Tehran, Iran 

Mehran Kardar 

Physics Department, Massachusetts Institute of Technology, Cambridge, MA 02139 

(Dated: February 1, 2008) 

The Casimir force between macroscopic bodies depends strongly on their shape and orientation. 
To study this geometry dependence in the case of two deformed metal plates, we use a path inte- 
gral quantization of the electromagnetic field which properly treats the many-body nature of the 
interaction, going beyond the commonly used pairwise summation (PWS) of van der Waals forces. 
For arbitrary deformations we provide an analytical result for the deformation induced change in 
Casimir energy, which is exact to second order in the deformation amplitude. For the specific case 
of sinusoidally corrugated plates, we calculate both the normal and the lateral Casimir forces. The 
deformation induced change in the Casimir interaction of a flat and a corrugated plate shows an 
interesting crossover as a function of the ratio of the mean plate distance H to the corrugation 
length A: For \ <^ H we find a slower decay ~ compared to the II~^ behavior predicted by 

PWS which we show to be valid only for A S> -ff. The amplitude of the lateral force between two 
corrugated plates which are out of registry is shown to have a maximum at an optimal wavelength of 
A « 2.5 H. With increasing H/\ > 0.3 the PWS approach becomes a progressively worse description 
of the lateral force due to many-body effects. These results may be of relevance for the design and 
operation of novel microelectromechanical systems (MEMS) and other nanoscale devices. 

PACS numbers: 03.70.-|-k, ll.lO.-z, 42.50.Ct, 12.20.-m 

I. INTRODUCTION 

More than five decades ago, Casimir predicted that the ground state energy of photons is alternated in the presence 
of two parallel perfectly conducting metal plates in such a way as to lead to an observable macroscopic force between 
them Q. The attractive force (per plate area A) has an universal amplitude, an energy scale set by the fundamental 
constant he, and decays with the distance H between the plates as 

F _ TT^ he 

A ^ ^ 240 ■ ^ ' 

This remarkable prediction of quantum electrodynamics has implications in many contexts ranging from surface 
physics Q, particle physics |^, to cosmology Q. Because of its fundamental nature, the Casimir effect has motivated 
extensive theoretical work, especially during the last decade. The pioneering result in Eq. (|l|) has been generalized to 
include important effects such as the finite conductivity and surface roughness of the plates, finite temperature, and 
even moving plates in the dynamic counterpart to the Casimir effect (see Refs. ||, 0, ||, ^, |ll[ for reviews). 

On the experimental front, the early attempts at observing the Casimir force, initiated by Sparnaay in 1958 
and later by Van Blokland and Overbeek in 1978 ]l3t , were not conclusive due to large experimental uncertainty. 
In recent years, however, there have been a number of precision measurements which set the modern stage in this 
field; starting in 1997 by Lamoreaux jl^ who used a torsion pendulum with an electromechanical feedback system 
to measure the Casimir force between a spherical surface (lens) and a flat plate. Mohideen et al. |jl^ measured 
the Casimir force between a sphere mounted on the tip of a flexible cantilever and a flat plate by an atomic force 
microscope. Chan et al. [l^ measured the Casimir force between a sphere and a flat plate in a microelectromechanical 
system (MEMS) using a micromachined torsional device. All these experiments conflrm the Casimir force formula in 
the range from 100 nm to several /im to a few per cent accuracy. In order to achieve this high precision, a careful 
analysis of the corrections due to the finite conductivity of the metal surfaces, roughness, and nonzero temperature. 
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is indispensable (Tl], [U, |T^, |T^, |2^, |2^. The above cited experiments in fact deviate from the fiat plate geometry 
corresponding to Eq.^ ) by using a plane-sphere configuration (thus avoiding the experimental difficulty of keeping 
two flat surfaces sufficiently parallel). The force for a spherical surface with (large) radius i? at a distance H of closest 
surface-to-surface approach from a flat plate can then be calculated from the Casimir potential £{H) for two flat 



plates by using the proximity force rule p^ ], see Eq. ( Pq ) below. Recently, however, G. Bressi et al. |23 measured 
the Casimir force between two parallel flat surfaces directly, confirming Eq. (0) to 15% accuracy. 

The Casimir force in Eq. (|l|) has analogies to the effective force between particles or plates immersed in a system 
close to its critical point, which arises due to the modification of thermal fluctuations of the bulk order parameter. 
This effect was originally predicted in 1978 by Fisher and de Gennes |Q for colloidal particles immersed in a binary 
liquid mixture near its critical de- mixing point, and observed experimentally for silica spheres immersed in a mixture 
of water and oil (2,6-lutidine) Related phenomena occur in liquid crystals |2^], microemulsions and for 

inclusions in fluctuating membranes pq] (see Ref. pC| | for a review). In recent years, the critical-point Casimir effect 
has attracted increasing theoretical [I^, ^if |3^, |3^ and experimental interest |3^, For He"' wetting 

films close to the superfluid phase transition, the theoretical predictions |^ for critical-point Casimir forces between 
parallel surfaces exhibiting Dirichlet boundary conditions have been confirmed quantitatively . 

Equation (^ for the electromagnetic Casimir force is valid in the ideal limit of perfectly conducting plates. In 
the more general context of the Lifshitz theory for dielectric bodies p6|] , this corresponds to an infinite dielectric 
constant e for all frequencies ui. For finite e — e{u!), this power law for the force is recovered for large distances 
H ^ c/ujq, where ojq is the smallest resonance (absorption) frequency of the dielectric (usually c/uiq « 10-100 nm). 
In this, so-called retarded, limit, the force is universal in the sense that it only depends on the electrostatic dielectric 
constant Sq = The opposite hmit of <C c/luq gives the unretarded van der Waals force F/A ^ . The 

interpretation of the Casimir force in terms of changes in zero point vacuum electromagnetic energy suggests it to be 
a strong function of geometry; probing the global shape of the boundary that confines the vacuum fluctuations 
Indeed, whereas the van der Waals force between electrically polarizable particles is always attractive, even the sign 
of the Casimir force is geometry dependent, and can be repulsive, e.g., for a thin spherical or cubic conducting shell 
psj ^ |4C| , |4l|] . Repulsive Casimir forces are expected also when magnetic properties of the boundaries are exploited, 
e.g., by using a perfectly conducting and an infinitely permeable plate 

Apart from the importance of these phenomena to basic science, the ongoing refinement of nanofabrication tech- 
nology in electronics and mechanics also provides new impetus for the understanding of such systems in view of 
applications in nanotechnology. On length scales of about 100 nm and below, the Casimir force becomes comparable 
or even dominant to other forces |Q , and thus must be taken into account in the design and operation of nanoscale 
devices Indeed, the experiment of Chan et al. |l^ demonstrates the possibility for novel actuation schemes in 

microelectromechanical systems (MEMS) based on the Casimir force. They also show that the Casimir force can be 
used to control dynamic properties of such systems, e.g., in a nonlinear micromechanical Casimir oscillator |^^. The 
above mentioned dependence of strength and sign of the Casimir force on geometry and material properties offers the 
opportunity to manipulate this interaction in a controlled way, e.g., by tailoring the shape of the interacting surfaces. 
On the other hand, movable elements of nanoscale devices may unwantedly stick together due to the strong attractive 



Casimir force, impeding their operation | |44| , 45 1. This so-called stiction could possibly be prevented by using shapes 



(e.g., suitable modulations) and materials of movable elements such that the Casimir force between them is repulsive 
(or at least reduced). 

Due to the importance of the Casimir force to basic and applied science, it is highly desirable to demonstrate its 
strong shape dependence in a set-up that clearly shows its distinction from the usual pairwise additive interactions 
p7| . In a previous Letter we pointed out that a promising route to this end is via modifications of the parallel plate 
geometry, since measurement of the putative repulsive Casimir interaction for a conducting sphere is experimentally 
difhcult. In searching for nontrivial boundary dependences, Roy and Mohideen | |49t examined the force between a 
sphere and a sinusoidally corrugated plate with amplitude a w 60 nm and wavelength A « 1.1 fiui. (This geometry was 
first suggested in Refs. .) Over the range of separations H w 0.1—0.9 /^m, the observed force showed clear deviations 
from the dependence expected on the basis of decomposing the Casimir force to a sum of pairwise contributions (in 
effect, an average over the variations in separations). This experimental result motivated our calculation of the exact 
Casimir force in the geometry depicted in Fig. |l|, without the assumption of pairwise additivity. Our analytic results 
[see Eq. (^) and Fig. |j hold to second order in a, and show that for fixed H the corrections due to corrugation 
strongly depend on A. In fact, for H/\ 3> 1 the correction is by a factor of H/\ larger than in the opposite limit of 
H /\ <ti 1 where the assumption of pairwise additivity is asymptotically correct. However, the experiments of Ref. 
are performed in the range of H/X « 0.1 — 0.8 where the corrections to pairwise additivity may not be significant 
enough to account for the observed deviations. In Ref. it has been suggested that these deviations are due to 
a lateral force that tends to preferentially position the spherical AFM tip on top of local maxima of the modulated 
surface (leading to a smaller separation and stronger force). Based on our results, we thus propose that the shape 
dependence of the Casimir force can be probed more clearly by going to modulations of shorter wavelength; a hard 
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but achievable goal. 

The originally predicted Casimir force between two flat metal plates (or between a flat and a deformed plate) is, 
for symmetry reasons, oriented normal to the surfaces of the plates. However, if both plates are deformed there is 
also a lateral Casimir force, as predicted in Refs. |^ and confirmed experimentally To date, the lateral force 
between two corrugated plates has been calculated explicitely only within the approximative approach of a pairwise 
summation of van der Waals forces, see, e.g., Here we calculate the lateral force exactly to second order in the 
corrugation amplitude for the geometry shown in Fig. ^, without referring to a pairwise summation scheme. As for 
the normal force, our results [cf. Eq. (Q) and Fig. ||] show that the lateral force strongly depends on the ratio of the 
corrugation length A and mean plate separation H . We find the pairwise summation to be a valid approximation only 
for sufficiently small values of H/\ < 0.3. However, the experiment of Ref. |52 is performed at a ratio H/X « 0.18 
where we do not expect significant deviations from the pairwise summation approximation. 

The outline of the paper is as follows: In the next section we set up the general path integral formulation for the 
Casimir energy of the electromagnetic gauge field. By separating into transversal electric and magnetic modes, we 
reduce the problem to two decoupled problems for scalar fields which differ only in their boundary conditions. In each 
case, we calculate the Casimir energy for general plate deformations perturbatively. In the third section, we give a 
brief summary of the pairwise summation approach, and the resulting Casimir interaction. Section four gives detailed 
results on both the normal and the lateral force between sinusoidally corrugated plates. We conclude by discussing 
the relevance of our results to experiments. Details of the calculations are left to the appendices. 



II. PATH INTEGRAL FORMULATION OF THE CASIMIR ENERGY 



We consider two perfectly conducting deformed plates Sa (a — 1,2) of mean separation H, which are infinitely 
extended along the plane spanned by y|| = (2/1,2/2)- Assuming static and uniaxial deformations without overhangs, 
their profiles are described by height functions ha{yi), with J dyi ha{yi) — 0. The Casimir energy at zero temperature 
corresponds to the difference of the ground state energies of the quantized electromagnetic (EM) field for plates at 
distance H and at ^ 00, respectively. To obtain this energy, we employ the path integral quantization method. 
For general deformations, it is necessary to consider the action 

^cm{A^} - - J I d'^XF^^F^^'^ , (2) 

where X denotes a point of 4D spacetime, F^j, = d^j^A^ — d^A^, and the four potential A^ is subject to the boundary 
condition that the tangential components of the electric field vanish on the surfaces. The redundant degrees of 
freedom due to the gauge invariance of the electromagnetic field can be eliminated by the Faddeev-Popov gauge fixing 
procedure ^ . 

However, for the uniaxial deformations under consideration here, we can develop a simpler quantization scheme, 
by a similar reasoning as used in the context of waveguides with constant cross-sectional shape [^^. In this case, the 
transverse magnetic (TM) waves and transverse electric (TE) waves (with respect to the translational invariant j/2 
direction) constitute a complete set of modes to describe an arbitrary electromagnetic field between the plates [ p4[ . 
For TM waves all field components are then uniquely given by a scalar function corresponding to the electric field 
along the invariant direction, 

$TM(i,2/i,y2,z) = E2{t,yi,y2,z), (3) 

with the Dirichlet boundary condition $tm|s„ = on each surface Sa- The TE waves are analogously described by 
the scalar function 

'^TE{t,yi,y2,z) = B2{t,yi,y2,z), (4) 

with the Neumann boundary condition c^h'I'teIsq = 0, where 9„ is the normal derivative of the surface Sa pointing 
into the space between the two plates. After a Wick rotation to the imaginary time variable X^ — ict, both fields 
$TM and 'I'te can be quantized using the Euclidean action 

^IJd^X (V$)' . (5) 

In the 4D Euclidean space, the plates are parameterized by ^i(y) = [y, /ii(j/i)] and ^2(7) = [y, + h2{yi)], where 
y = (2/0,2/1,2/2) = (2/0, y||), and 2/0 = ict. 
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In order to obtain the ground state energy from this quantization scheme, we now consider the partition functions 
and Zjvf for the scalar field Euclidean action both with Dirichlet (D) and Neumann (N) boundary conditions 
at the surfaces. Following Refs. ^ (cf. also Ref. Q), we implement the boundary conditions on Sa using delta 
fimctions, leading to the partition functions 

1 r ^ 

= ^ n^^I[I[^M^c.)]eM-SEm/h), (6a) 

1 r ^ 

2n = — V^l[l[S[dn-fiX^)]expi~SEm/h), (6b) 

^ ^ ^ TV 



a=l Xa, 



where Zq is the partition function of the space without plates. Note that in Eqs. ( paj ) and (6b) no gauge fixing 
procedure is needed, since the relevant degrees of freedom are expressed in terms of the fields E2 and B2 itself rather 
than the vector potential in Eq. (|^). The details of the calculation of the partition functions are left to Appendix 
0, and yield 

ln2^D = -^trlnEo, ln2^N = -^trlnEN . (7) 
The kernels Ed and En are given by 

[rD]a/3(y,y') - [5.(2/1)]^/^ G[X„(y) - X^(y')] [5My'i)]^/\ (8a) 

[T^Uiyy] = K(2/i)]'/'5„„(,,)5",fo;)G[^a(y)-^/3(y')][5My;)]'/\ (8b) 

where ga{yi) = 1 + [h'a{yi)Y is the determinant of the induced metric, and na{yi) — (— l)".gQ ^^^(j/i)[/i„(?/i), 0, —1] is 
the normal vector to the surface Sa, while 

is the free correlation function (Gaussian propagator) corresponding to the Euclidean action S'e in Eq. (|). In the 
following, we denote position vectors in the 4D Euclidean space by r = (y, z). We can then extract the Casimir energy 
£ per unit area as 



£{H) = E{H) - lim E{H) , (10) 

H^OQ 



with 



he 

E{H)^^—[\nZ^+\nZ^\, (11) 

where A is the surface area of the plates, and L denotes the overall Euclidean length in time direction. 

The above equations provide an exact result which yields the Casimir energy for arbitrary static uniaxial deforma- 
tions. These equations can be used to evaluate the Casimir force by a recently developed numerical approach [E7[. 
However, in order to obtain a closed analytical expression for the partition functions in terms of the kernels in Eqs. ^), 
here we resort to a perturbative expansion with respect to the height profiles ha{yi)- In fact, for this expansion to 
be valid, we have to assume that the amplitude of the deformations sets the smallest (geometric) length scale of the 
system. In what follows, perturbation theory is carried out to second order in ha{yi), separately for the two types 
of boundary conditions. However, within second order in the height profile, the result is exact in the sense that it 
correctly takes into account the many-body nature of the Casimir interaction. 



A. Dirichlet boundary conditions 

Following Refs. we expand InZo in a series IuZdIo + Ih-^dIi + lnZD|2 + ■ . ., where the subscript indicates 

the corresponding order in ha. The lowest order result is 
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corresponding to two flat plates. The first order result Ih-ZdIi vanishes since we assume, without loss of generality, 
that the mean deformations are zero, J dyihaiui) = 0. The complete second order contribution is given by 

in^Db = -IJd'y miyi)? + [h'Mr] / (03 X 1 (13) 

+ 1 J d'y J dV/^D(|y-y'|) {^-llhiiyi) - hM)? ~ l[h2{yi) ~ hMf^ 
-IJ'^'yJ d'y'QDi\y-y'\) MyMy',) + hMh,{y[)] . 

In the first term, which is further discussed below, h'^ = V/iq,. The kernels appearing above are given by ||59|| 

KB{y) = Fi{y)d^G{y,0) + Fiiy)F5{y) + F^{y)\ (14a) 
QD(y) = FMdlG{y,H) + F^F^iy) + F2{yf , (14b) 



with the set of functions 



- J^--W^y 

where y = |y|, p = |p|, G{p,z) ~ ^e^^'^l is the partially Fourier transformed free propagator of Eq. (||), and 

N{p, H) = [G(p^ Q)Y — [G(p, H)]^. The functions in Eqs. ( p^ ) can be calculated explicitly (see Appendix [B l| ), leading 
to the simple result 

ifD(y) = 2i^D,oo(?/)+i^D,rog(y), (16a) 

Qj3(y) ^ 1 (16b) 

^"^^^ 128 H^y^ cosh\s) ^ ' 

with 

1 1 TT^ 1 cosh^(s) 

ifD,oo(y) = ^, = + ^> (17) 

where s — 7ry/(2iJ). The kernel i^D has two contributions of different origin. In the limit H — > oo, corresponding to 
two decoupled surfaces, one has K^) — > 2_R'd,oo, while Qoiy) vanishes. Thus the part /Cd^oo describes a single surface. 
The two (/^-independent) single surface contributions have to be subtracted from the total kernel Kd in order to 
obtain the regularized kernel i^o.rcg which has to be used in the calculation of the Casimir energy in Eq. (UQ). For 
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finite iJ, the kernel K^{y) actually has contributions from both outside and inside the cavity, whereas i^D.oo comes 
from outside and the second term of ifD.rcg from inside. The kernel Quiv) has only contributions from inside the 
cavity. 

It is instructive to discuss the meaning of the contributions to lnZD|2 in Eq- ([l3|). The terms in the first row are 
iJ-independent and formally divergent. They do not contribute to the Casimir force between the surfaces but yield a 
quantum electrodynamical increase of the surface tension of the individual surfaces after introducing a suitable short- 
distance cutoff [ pO[ . The necessity for a cutoff stems from our continuum approach which breaks down on microscopic 
length scales. The remaining terms in Eq. jl^ ) all contribute to the Casimir force (with Kd replaced by ifD.rcg)- The 
local contributions in the second row are half (due to TM modes only) of the individual surface (non-mixed) terms 
which follow in second order of perturbation theory in h^, from the pairwise summation approach, cf. the second term 
in Eq. (p^). The third row in Eq. ( p^ describes non-local individual surface contributions which are missing in the 
pairwise summation approach. Finally, the last row accounts for contributions due to the interference between the 
two surface profiles. Obviously, it has a more complicated form than the corresponding last term in the approximative 
pairwise summation result in Eq. (p6|). 



B. Neumann boundary conditions 

Expanding IuZn in a series with respect to ha as before, the lowest order result is the same as for the Dirichlet 

case, 

, ^ , AL tt'^ , 

and the first order result InZivjIi again vanishes. The complete second order result assumes a similar form as for 
Dirichlet boundary conditions. We find 



In^Nb = +ljd'y {[h[iyi)? + [h'Mr} / X 1 (19) 



240 



+ \ J d'y J d'y'K^{\yo-y'o\Ayil-y\\\)[-l[hi{yi)~h^iy[)r~^[hM~hMr^ 
-i / dhj j dVQN(l2/o-2/oMy||-y||l) [KiyMy'i) + hMh^{y[)] . 



The kernels for Neumann boundary conditions assume a more complicated form since the normal derivative breaks 
the equivalence of space and time directions. The result reads 

^N(lyoUyiil) = Ti{y)dlg{y,0) + T,{y)H{y) + J^^{yf (20a) 
+ d,d,[Ti{y)d,d,G{y,Q) + J'liy) d,djT^{y) + d^J'^iy) djT^{y)] 
+ 2d,[J^^{y)d,g{y,Q) -f Ti{y) d.J'^^iy) + T^{y) d,T^{y)] , 

QN(lyoUyiil) = ^i{y)dlg{y.H) + Ti{y)T(,{y) + T^{yf (20b) 
+ didj[J^A{y) didjG{y,H) + J^i{y)didjT8{y) d,Tw{y) d^Tw^y)] 
+ 2^dTM^^g{y.H) + Ti{y)d,J'M + ^2(y) a,.Fio(y)] , 

where summation over i = 1,2 and j = 1,2 is understood. Note that di and dj act on the spatial components yy 
of y only. This is the reason why the rotational symmetry within the 3 dimensional y space is broken for Neumann 
boundary conditions. For translationally invariant profiles ha{yi) only terms with i — j — \ contribute as can be 
seen by integration by parts. However, as in the case of Dirichlet boundary conditions, the above result is valid for 
any /ict(y) but then can no longer be interpreted as the contribution from TE modes to the electrodynamic Casimir 
energy. 

Here, we have introduced 

9{y,^) - / (|^e^P '^g(p,z), (21) 
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with g{p, z) = d^G{p, z) = |e ^'^1, and the functions J^j{y) defined as 



My) = 
My) = 
My) = 
My) = 
My) = 
My) = 
My) = 
My) 

M{y) 

Miy) 

^i2{y) 



(27r)3 ' 


0) 

viP,H) 


dzg{p, H) , 


^'P e^p y 
(27r)3 


9{P,H) 
i]{p, H) 


dzg{p, H) , 


(27r)3 


9{P, H) 
i]{p, H) 


1 


d^P rp.y 

(2^)3 


5(P,0) 




^"^^ „«p y 
(27r)3 ^ 


9{P,H) 
i]{p, H) 


[d.9{p,H)f 


'^'P e^P y 
(27r)3 


9{P,0) 
i]{p, H) 


[d.G{p,H)f 


d'p ,p.y 

(2^)3 


9{P,H) 
i]{p, H) 


[dzG{P.H)f 



d^p gjp.y gjp, H) 

(27r)3 ij{p,H) 



d^P ^p.y g(p,o) 



(27r)3 



7?(p,i7) 



d.G{p,H), 



dzG{p,H), 



7 (27r)'^ r]{p,H) 



(22b) 
(22c) 
(22d) 
(22e) 
(22f) 
(22g) 
(22h) 
(22i) 
(22j) 
(22k) 
(221) 



with r]{p, H) = [g{p, 0)]^ — [(7(p, H)]'^. The expUcit form of these functions can be found in Appendix [B 2| . 

The result in Eq. (19) has the same type of contributions as discussed for the Dirichlet case. Both Dirichlet and 
Neumann case include 'surface tension' contributions, but with opposite signs, and identical local terms (second row 
in Eq. (|l9|)). Since these local terms are the only (non- mixed) contributions obtained by the pairwise summation 
approach, the latter does not distinguish between the two types of boundary conditions. The main results of our 
general analysis of surface deformations are contained in Eqs. (13), (19). In section IV we apply these results to the 
important case of modulated plates. 



III. PAIRWISE SUMMATION APPROXIMATION 



The path integral approach may be compared with the commonly used approximative method of pairwise summation 
(PWS). In the latter approach, the Casimir energy £(iJ) for two arbitrary shaped bodies of mean distance H is obtained 
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by the pairwise summation of a two-body potential U{r). In terms of the deformation fields this leads to 



£{H) = - I dVll / d'y'n / dz I dz' U 

J H+h2{yi) 



hi(yi) 



— oo 



1/2' 



(23) 



In general, these integrals need to be computed numerically. However, there are the following simplifications. If one 
of the plates is flat, e.g., hi{yi) = 0, the integrals can be performed explicitly, leading to the simple result 

£iH) = ^l d'y\\£o[H + hM], (24) 

where Sq is here the energy of two flat plates at distance H, calculated from the same pair potential U{r). Thus, 
in this particular case the pairwise summation approximation is equivalent to a geometrical average of the flat plate 
energy with locally varying plate distance over the plate area. 

For two deformed plates, the integrals in Eq. (|2^) in general can only be performed perturbatively in the height 
profile. To do so, we follow the usual PWS approximation and assume a 'renormalized' retarded van der Waals 
potential ||ll|, psf , 

U{r) = -^-^r-\ (25) 

The 'renormalization factor' of the pair potential is chosen here such that in the limit of two flat plates the exact 
Casimir result £q, cf. Eq. (|3^), is recovered. To second order in ha one obtains now 

(26) 



^ ' 720 m 120 H'^ A J •'^ll L ivyi; "T 2vyiyj ^ y y\\ j y 



(i?2+y'2)7/2 



For simple types of plate modulations, the integrals over the deformation fields can be calculated easily. This will 
allow us a direct assessment of the validity range of the pairwise summation approach by comparing Eq. ( p6[ ) to the 
predictions of the path integral technique. 

IV. MODULATED PLATES 

We now apply the results of Sec. || to static uniaxial modulations of two parallel plates. In the first part we focus 
on the normal Casimir force between a flat and a corrugated plate. This force per unit area is defined as 

F.^-H, (27) 

in terms of the Casimir energy in Eq. ([l0|). However, in most of the experiments, the flat plate is replaced by a 
spherical lens with large radius R ^ H. In the latter case the normal force can be obtained by using the Derjaguin 



approximation (DA) (or proximity force rule) |22 , leading to 

^^DA,n = 2TrR£. (28) 

Therefore, in the context of the normal Casimir force, we just calculate £ explicitly. 

In the second part of this section we generalize our results for two modulated plates with equal modulation length 
but with a phase shift between them. Due to the broken translational symmetry, there is now also a lateral force 
between the two plates which arises solely from the cross-terms /ii/i2 in Eqs. (jl^), ([l9|). If we denote the shift 
between the two corrugations by the length b, the lateral force is obtained from 

A. Normal force 

As a prototype of a corrugated surface, and to make contact with recent experiments by Roy and Mohideen p9[ , 
we consider a sinusoidally modulated plate along the yi direction, with amplitude a, wavelength A, and mean distance 
H from the flat plate (see Fig. |l|), i.e., 

/ii(yi) = acos(27ryi/A), and /i2(yi) = . (30) 
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FIG. 1: The set-up used for calculating the Casimir energy between a flat and a corrugated plate at mean separation H. 



For this particular deformation [Q, only the Fourier mode of wavelength A in the kernels K-oly) and i^TNdyol, |yi||) 
is probed. Thus, the calculation of InZo and IuZn via Eqs. (13), (|l9|) reduces to Fourier transforming the kernels 
with respect to yi. The corresponding expression for £ in Eq. (IC) can be written as 



f + ^ci ; 



with 



TT^ he 

720 



£cf = — 



Gtm 



Gte 



(31) 



(32) 



where Eq is the energy of two flat plates and the index cf of £cf stands for corrugated-flat geometry. The notation 
0{a^) indicates that third and higher powers of a/H and a/A are not considered here. The corrugation induced 
contributions to £ci from TM and TE modes at second order in a are governed by the functions 



GTuix) = — -t- go{x) 



Gte (a;) 



480 



gi{x) +xg2ix) +x^g3{x) 



(33a) 
(33b) 



The first term 7r^/480 in both equations corresponds to the local contributions, cf. the second row in Eqs. (p^), 
(|l9|). Since these are the only terms which are obtained within a pairwise summation approach, the functions gm{x) 
represent non-trivial corrections which are neglected in the pairwise summation scheme. These functions can be 
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calculated from the kernels Ku^y) and i^Ndj/ol, |y|||) by Fourier transformation, leading to the expressions 



TT 

480 
480 

2 



TT 

480 ■ 
7r3 



TT 

128 



°° , sin(4a;s) 
as ■ 



64 



1440 

sin(4a;s) 



4xs 



13 

4^6 



4xs 



5 

3^ 



cosh^(s) 



15s' sinh^(s) 



(34a) 
(34b) 



4 5 cosh(s) 
455^ 2s^ sinh(s) 



3 cosh(s) 
2^ sinh^(s) 



+ 



64 



1 cosh'^(s) 1 cosh(s) 1 cosh^(s) 5 cosh^(s) 1 1 
2^ sinh^(s) ^ " ^ " ^ ' " " 

ds cos(4xs) 



s sinh^(s) 2 sinh''(s) 4 sinh''(s) sinh^(s) 



52(2;) 



TT 

64 



TT 

32 



ds sin(4a;s) 



2 
45^ 



45s2 
5 

7E + 



2 

3^ 



11 1 cosh(s) 1 cosh(s) 

s4 sinh2(s) ^ sinh(s) ^ sinh^(s) 



1 cosh(s) 



s'^ sinh (s) 



ds sin(4xs) 



s'^ sinh(s) 

1 cosh(s) 1 cosh'^(s) 
sinh^(s) s sinh'*(s) 



(34c) 



4 
3s 



2 cosh(s) 
sinh(s) 



1 cosh^(s) 
s sinh^(s) 



(34d) 



Before giving the explicit forms of these functions, let us consider two limiting cases. For a; — s- 0, or A/iJ — s- 00, we have 
go{x), gi{x) 0, and 52(2;), gsix) converge to finite numbers, thus leaving in Eqs. (|3^ ) only the local contributions 
from the pairwise summation approach. In the opposite limit x —^ 00, or A/if — *■ 0, the integrals in Eqs. ( |3^ ) decay 
to zero, leading to 



50 (a;) 



TT 

480 



TTX ■ 



1 



5tt 1 
126x 



TT 

480 



1 



TT 1 

18x 



0{e- 



(35) 



and both 172(2^) and 53(2;) are ©(e"**^^). From this result it is obvious that for \/H ^ in Eqs. (p3|), the terms 
7r^/480 from the pairwise summation approach are exactly canceled by corresponding terms of opposite sign in the 
non-trivial corrections described by go{x) and 51(2;). The most relevant contributions in this limit arc now provided 
by the first term in 50(2^) and 51(2;) in Eq. (^5|), leading to the novel scaling behavior Gxm/te(2;) ~ 2; + 0{l/x). 

The integrals in Eqs. ( ^^ can be carried out for A > 0, or equivalently a; > 0, by closing the integration contour 
via a semi-circle at infinity in the upper half of the complex s plane, using the residue theorem pl| |. The resulting 
sum of an infinite series of residues can be expressed in terms of the polylogarithm function 162] 



leading to, with u = exp(— 47ra;) 
Gtm{x) 



480 



+ 



647r 



30 

2Li5(w) + 



ln(l - u) 
1 



2567r3a 



Li„(z) 



E 



1920a; 
Li6(u) 



Li2(l-u) 

6 



TT 

945 



TTX^ 2 

Li2(u) H Li3(u) H Li4(M) 



(36) 



(37) 



Gte(2:) 



TT^'X TT'-X^ , ,^ , T- • /, N ^2; 

ln(l - u) H Li2(l - u) (1 

1440 30 ^ ' 1920a; ^ ' 48 ^ 



1 



2567r3a 



2a;^) Li2(u)+ 7^-777 ]l-h{u) 



48 64 



7 TT*^ 
-Uq{u) - TT'^Ui{u) + — 



(38) 



It should be noted that the appearance of the polylogarithm function in quantum electrodynamics is also known from 
the fine structure constant dependent corrections to the gyromagnetic ratio of the electron |63]. 
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Figure |^ displays separately the contributions from Gtm and Gte to the corrugation induced correction Set to the 
Casimir energy. While Gtm{H / X) is a monotonically increasing function of H/X, G^EiH / X) displays a minimum for 
H/XkO.3. 




LI OJ I'l-I IJ.M its J 



iifj. 

FIG. 2: Rescaled correction Ed to the Casimir energy due to the corrugation as given by the terms in square brackets of 
Eq. ( ^l| ) (upper curve). The lower curves show the separate contributions from TM and TE modes. The rescaling of S^i is 
chosen such that the corresponding prediction of the pairwise summation (PWS) approximation [second term of Eq. (^6|)] is a 
constant (dashed lines). 




15 10 15 

H/a 

FIG. 3: Rescaled Casimir energy to second order in a/H and a/A as given by Eq. ( |3l| ) for two fixed values of a/A, shown as 
dashed curves. The rescaling is chosen such that the Casimir energy of two flat plates becomes one (horizontal line). The curve 
for a/X — 0.05 corresponds to the parameters used in the experiment of Ref. [jlSJ], where H/a varies between approximately 3 
and 17. Note that the dashed lines are accurate predictions for the true Casimir energy in the limit a <^ H,X only, as indicated 
by the question mark (see text and the discussion in Sec.^). For comparison, the result of the pairwise summation (PWS) 
approximation [cf. Eq. (|2^)] is shown. It agrees with the perturbative result in the limit A = oo only. In the opposite limit. 
A/a — > 0, the energy can presumably be estimated from a reduced distance argument ^ (see Eq. (pl[)). 

Examining the limiting behaviors of Eq. ( pl| ) is instructive. In the limit X ^ H, the functions Gtm and Gte 
approach constant values, and the Casimir energy takes the A-independent form 

Note that only in this case both wave types provide the same contribution to the total energy and the result agrees 
with the pairwise summation approximation (see Fig. 0). In the opposite limit of A ^ iJ, as demonstrated above, both 
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Gtm and Gte grow linearly in H/ A. Therefore, in this limit the correction to the Casimir energy decays according 
to a slower power law in , as 

^ = -|^^(' + ^4) w 

with an amplitude proportional to 1/A. Note that this behavior is completely missed by the pairwise summation 
approach which always yields a A independent Casimir energy in the presence of modulations on one plate. 

In the limit A ^ a, {H — a) we expect that the factor multiplying a/ H in Eq. ( ^0| ) saturates at a number of order 
unity. This result can be justified by noting that the most relevant contributions to the force come from modes of 
wavelength of order H. The corrugation also affects modes of wavelength of order A, but these modes contribute to 
the single plate energy only. Thus, in the extreme limit A ^ a, {H — a), one has a clear separation of the length scales 
H and A, and the modes "see" flat plates at the reduced separation H — a, (a > 0) Esl. More recently, an exact 
approach for calculating the Casimir force has been developed which confirms the above argument and yields for the 
case of TE modes (Dirichlet boundary conditions at both plates) the exact result , 

This leads to a correction of the order a/H (with prefactor 3) after expansion in a. 

The above behavior of the correction £cf for small and large H/X clarifies the limits of validity of previous results in 
the literature. The upper dashed line in Fig. |^ corresponds to the PWS approximation (see Refs. ||ll|, Q and Sec. Ill) 



It is evident that this approximation is accurate only for H/X (which in this limit is equivalent to the Derjaguin 
method to any order in the amplitude a Already for H/X of order unity, the PWS approximation breaks down. 



The opposite limit, H/X oo, corroborates the result reported in Rcf. |17 , which is larger than the former by a 
factor of H/X ^ 1. However, in experiments with lateral distortions A of the order of H, none of the above limiting 
cases is realized, which makes the present, more complete analysis necessary. 

The use of a spherical tip, of large radius R, in experiments causes some differences from the flat plate geometry 
used in our calculations. First, the positioning of the tip relative to the modulations is important when H and A 
are comparable, but becomes insignificant in the proposed limit of A <C iJ, R- Secondly, as long as R ^ H,X the 
curvature of the tip does not lead to nontrivial corrections, and the force can be related to the energy per surface 
area £ in Eq. (^) by the proximity force rule F = 2itR£ in Eq. (|2^). These formulas thus provide a specific recipe 
for evaluating the nontrivial shape dependences of the Casimir force in the experimental set-up. 

The net Casimir energy £ is shown in Fig. ^ for two representative values of a/ A, including the parameters used in 
the experiment of Ref. Note that the corrugation induced correction leads to a larger energy £, and hence the 
corresponding (attractive) force F — 2t:R£ is enhanced, at least to second order in a/H and a/ A, which becomes exact 
in the limit a <^ H,X. This trend suggests, in particular, that in the set-up of Fig. ^ the force is always attractive, 
although definite statements for values of a/H and a/X of order one can only be made by using non-perturbative 
methods, as indicated by the question mark in Fig. || (see Ref. |57 and the discussion in Sec.^). However, in the 



opposite limit, for which the tips of the modulations of the lower plate in Fig.|l| almost touch the upper (flat) plate, 
i.e., H — a <^ A^/a, a, the energy can be calculated exactly by using the Derjaguin approximation for the individual 
tips of the modulations; this leads to 

which corresponds to the result in Eq. (^ ) after taking H — a ^ 0. The above result implies that at least for the 
particular case of an uniaxial sinusoidal corrugation the corresponding force F = 2ttR£ is attractive when the surfaces 
almost touch. 

B. Lateral force 

As a natural generalization of the geometry of the previous section, we study the Casimir interaction between two 
sinusoidally corrugated plates. For direct correspondence to recent experiments [Q, we consider the specific profiles 

hi{yi) — acos{2nyi/X) , and /i2(yi) = acos (27r(?/i -I- 6)/A) , (43) 

which arc shifted relative to each other by the length b (see Fig. 
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By inspection of the deformation dependent contributions to the Casimir energy in Eqs. (y_3|) and (|19|), one obtains 
for the total Casimir energy £ of the corrugated-corrugated geometry the relation 



(44) 



with £q and £cf given in Eq. (p^), and the corrugatio n-cor rugat ion i nteraction energy £cc which can be calculated in 
terms of the kernels Qriiy) and QN(|yo|, |y|||) in Eqs. (14b) and (20b). Besides oscillating contributions to the normal 
Casimir force Fn{b) from £cc{b), a lateral force 



d£cc 
db 



(45) 



is induced by the corrugation-corrugation interaction. This lateral force is much better suited for experimental tests 
of the influence of deformations, since there is no need for subtracting a larger baseline force (the contribution of flat 
plates) as in the case of the normal force. 

The calculation of the interaction energy fee again reduces for sinusoidally corrugated plates to Fourier transforming 
the kernels Q-oiy) and (5n(|2/o|j |y||l)- Separating the contributions from TM and TE modes, we find 



hca^ 



with 



Jtm{x) = ja{x), 

Jte{x) ^ - x'^ j2{x) + x"^ jsix) 



and the functions jm{x) are given by 

joix) = 

jiix) 

hix) 

jsix) 



32 



ds 



TT 

32 



= — / ds 



T 

TT^ 

2 



ds 



ds 



sin(4xs) sinh^(s) 

4xs cosh^(s) 

sin(4xs) sinh^(s) 

4:xs cosh^(s) 

sin(4xs) sinh^(s) 

4xs cosh'^(s) 

sin(4xs) sinh^(s) 



sinh (s) 

2 ^ ^ 



4xs cosh (s) 



(46) 

(47a) 
(47b) 

(48a) 
(48b) 
(48c) 
(48d) 



Before giving explicit forms for these integrals, it is instructive to again consider their extreme limits. For x = H/X — s- 
0, we find that both functions Jtm{x) and Jte(2;) tend to 7r^/240. In the opposite limit x = H/X —> 00 both functions 
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■0.1 a' 




FIG. 5: Rescaled lateral Casimir force amplitude \Fi as obtained from the path integral approach [Eq. ( |50[ )] (solid curve), and 
from the pairwise summation approach (PWS) [Eq. (5l|)] (dashed line). The results hold to second order in a (cf. the discussion 
in Sec.0). 



decay exponentially fast to zero so that the lateral force vanishes in this limit. In order to get the behavior in between 
these extremes, we have to calculate the integrals in Eqs. (^). Using the residue theorem, we finally obtain after 
summing over an infinite series of residues by using the Lerch transcendent [|63 



$(z, s, a) 



y " 



(49) 



the results (with u — cxp(— 47r.T)), 



^TM(a;) = 



120 



[l6x^ — l) arctanh(-y?I) + Vu 



IT 

12 



1 

80x 



16^*^"'^' 2) 



1 



— ^$(?/,5,-) + — ^$(?/,6, 



12 '^("'^'2) 



(50a) 



^TE(a;) = (iSa;"* - 1) arctanh(\A«) + 



17 17 1 
^ I X ) $(m,4, -) H $(u,5, -) H 5- $(ii,6, -) 



(50b) 



This result can be compared to the pairwise summation approach by considering the last term in Eg. (|2q). For the 
surface profiles considered here (cf. Eq. (^)), this term provides an interaction energy given by Eq. ( ^6| ) with the sum 
Ji:yi(x) + Jte{x) replaced by the function 

■hw^x) = — {Att^x'^ + 6t:x + 3) V^. (51) 
360 

The two results agree for x = H/X 0, since Jtm(O) + Jte(O) = 7r^/120 = >/pws(0). At the other extreme of 
A ^ iJ, both Jtm{x) + Jte{x) and Jpws{x) decay exponentially fast but with different H/X dependent coefficients. 
In particular, for large x = H/X, we get to leading order 

47^2 

Jtm{x) + Jte[x) = — {x^ + 0[x''))^/^ (x^oo), (52) 
to 

in contrast to the ~ x"^ behavior in Eq. (pT|). 

Since both Jtm {x) + Jte {x) and Jpws {x) are positive for all values of x, the equilibrium position of two modulated 
surfaces is predicted at 6 = A/2 in both approaches. This corresponds to aligning the maxima and minima of the 
two corrugations (cf. Fig. 0). The amplitude of the lateral force per unit area [Eq. (^5|)] as obtained from both 
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approaches is plotted in Fig. g. Interestingly, for fixed H there is an optimal modulation length A « 2.5H at which 
the lateral force is largest. Our result shows that the pairwise summation approximation is not justified beyond 
H/ A w 0.3. For H/ A of order one, the pairwise summation approach has already a relative error of about 150%. With 
increasing H/ A this error grows monotonically. 



V. DISCUSSION AND OUTLOOK 

We calculated normal and lateral Casimir forces between perfectly conducting modulated plates (Figs.|l| and ^) 
by means of the path integral quantization method (see Refs. |5^, and Sec.||). Based on the resulting exact 
expressions for the Casimir energy (Eqs. (|^) - (0)), we performed a perturbative calculation to second order in the 
deformation parameter a to obtain the results outlined in Sec. IV and shown in Figs. || and || (normal force) and Fig.|^ 
(lateral force). These results are thus exact to second order in a, and correctly take into account the many-body 
nature of the Casimir interaction, going beyond the commonly used pairwise summation (PWS) of van der Waals 
forces pi] , |5^ . Our results show significant deviations from PWS to second order in a. However, for finite values of 
a/ A and a/H, there will be corrections at higher orders in a/ A and a/H; in the present experiments, the sensitive 
range of a/H is of the order of 0.2 while we suggest values for a/ A of order one to probe the nontrivial shape 

dependence of the Casimir force (see below). At present, it is not clear how relevant the perturbative results to second 
(or any finite) order in a are for this range of the parameters a/ A and a/H (cf. the related discussion in Sec. |IV A[ ). 

To make further connection between our findings and the experimental situations, corrections due to the finite 
conductivity of the plates, finite temperature, and surface roughness should be taken into account as well ||l^, |l^, 
|l8| , [l9| po| , pl[| . These corrections introduce additional length scales into the problem, which are in turn the plasma 
wavelength Ap of the plates (e.g., Ap ~ 100 nm for aluminum [Q), the thermal wavelength At = hc/ksT (« 8 /xm at 
300° K), and the transverse correlation length ^ of the roughness (usually ^ « 300 nm jl^l). While the effect of finite 
conductivity becomes important for _ff < Ap, finite temperature starts to affect the result when H > Xt- The plasma 
and thermal wavelengths thus provide lower and upper bounds for H, respectively, such that our results for perfectly 
conducting plates at zero temperature are valid for X, H ^ Xp, and H <^ Xt- 

The importance of stochastic surface roughness can be deduced from our calculations. The relative corrections 
£ci/£ to the Casimir energy due to roughness of amplitude a and transverse correlation length ^ should be of the 
form (to second order in a) a^ /H^ for ^ ^ H, and a'^/{^H) for ^ ^ H. The latter behavior is in accordance with 
Ref. [y_8|. The experimental case corresponds to neither extreme, making a more complete analysis necessary. 

Given the above mentioned limitations, as well as the technical difficulties in achieving the desirable geometries in 
experiments, it seems that it is difficult to conclusively establish the nontrivial boundary dependence of the Casimir 
force. Figure ^ shows that in the experiment of Roy and Mohideen |Q, the lengths a and A are such that the 
nontrivial dependence of the Casimir force on the boundary shape is rather weak within the monitored range of H/a, 
and a pairwise summation of two body forces is a possibly adequate approximation. Our results suggest that a set-up 
with A of the order of a is better suited for observing the nontrivial geometry dependence predicted above. 

In general, one expects that as long as the nontrivial features in the geometry of the plates appear only as small 
perturbations to the trivial flat-plate geometry, the corresponding many-body effects of these features will be hard to 
measure. Considerably larger effects could result, however, in patterned surfaces with geometrical features that come 
close together across various parts of the surfaces. In such circumstances a non-perturbative calculation of the forces 
becomes necessary. Indeed it is most desirable to find robust numerical schemes (possibly along the lines of Ref. ) 
that can also incorporate the finite conductivity effects and surface roughness typical for experimental set-ups. 

Finally, we note that in the set-up of Fig. |^, nontrivial shape dependencies appear as corrections to a larger Casimir 
force. For the purpose of experimental tests, it is much more desirable to devise set-ups which directly probe differences, 
without the need for subtracting a larger baseline force. For example, in an atomic force experiment, simultaneous 
scanning of a flat and corrugated substrate would be desirable; while in the torsion pendulum experiment, one can 
imagine suspending a spherical lens equidistantly from two plates, one of which is corrugated. 
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APPENDIX A: PATH INTEGRAL FORMULATION FOR PARTITION FUNCTIONS 

We consider N manifolds (objects) with a = 1, . . . ,N. Each point on the manifold fla is represented by a vector 
^a(y) = (^a(y); M — ^ D dimensional manifold fla embedded in d dimensional space is parameterized by 

y = {yi, ■ ■ ■,yD)- 



1. Dirichlet boundary conditions 

The Dirichlet boundary condition $ = on the manifolds, can be enforced by the functional Ylx ^i^i-^a)] in 



Eq. (|6a|), which can be expressed in terms of auxiliary fields ^'q(Xq) as [47, pq 



[]<5[$(X,)] = y"l?*„(X„)exp I j dX^^^{X^)<^{X. 



(Al) 



The Gaussian integration over $ in Eq. (pa) can then be performed, resulting in 



N 



2d = / n^*"(^" 



,-Seff{*} 



The effective action S'off is given by 



(A2) 



(A3) 



a/3. 



where G{r,r') is the two-poi nt c orrelation function in unbounded bulk space. The functional integration over the 
curved manifolds ila in Eq. (A2) is facilitated by introducing the new fields ipa{y) = \ ['a[^a(y)1- However, this 
transformation requires some care regarding the integration measure dX^ in Eq. (A3), as well as the functional 



measure / ^^^(Xq) in Eq. (jAJ). The result is ||56|, |6J 

n2?*a(Xa)e-^-{*> = 



yje 



(A4) 



where the field 0a (y) = [5Q(y)]^^'*'0a(y) is given for each manifold Qa in terms of the determinant ga{y) of its induced 
metric 



5a,»j(y) = > -t: — ^ — ■ 
dyi dyj 

The new effective action ScS is then given by 

a/3 

where 

ra/3(y,y') = [5a(y)]^/*G[X4y),X^(y')][.9My')]'/' ■ 



(A5) 



(A6) 



(A7) 



The functional measure / l?0a (y) on the right hand side of Eq. ( A4 ) is the one conventionally used on a flat manifold 
(the local coordinate system). T he c orresponding Gaussian integrations can thus be performed, resulting in Eqs. (0) 
and (js^), with Fd = F from Eq. (A7). Note that in the present formulation, the trace and products of F are carried 
out by integrating y over the flat manifold of the local co ord inate system. Any dependence of Zd on the metric 
5a,ij(y) is contained explicitly in the definition of F in Eq. ( |A7| ). 
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2. Neumann boundary conditions 



For the Neumann boundary condition c}„$ — on the manifolds, the boundary condition enforcing functional 
Wx ^\'^n'^^^a)\ Eq. (6b) can again be expressed in terms of the auxiliary fields "^ai/^a) as 



[],5[a„$(X„)] EE yi)*„(X„)exp iy"dX„*„(X„)a„$(X, 



(A8) 



(A9) 



where the second line follows from an integration by parts. The Gaussian integration over $ in Eq. (6b) can then be 
performed, resulting in 



N 



with the effective action 



(AlO) 



(All) 



Sesm = ]-Y, I dXo, I dXp [dn^^(X^)]G{X^,Xp)[d,,^p{Xp)] 

= \ Y1 j dX^ j dXp-^^{X^)[dn^dn,G{X^,Xp)]^p{Xii) 
Calculations along similar lines as in the previous paragraph then lead to Eqs. (^ and (8b). 



APPENDIX B: CALCULATION OF THE KERNELS 

1. Dirichlet boundary conditions 

The kernels for Dirichlet boundary conditions were defined in Eqs. ( p^ in terms of the functions in Eqs. (p^. The 
explicit form of these functions is given by 



Fiiv) 



F2{y) 



1 TT cosh(s) 



8il3y sinh3(s) ' 



TT sinh(s) 



IQH^y cosh3(s) 



FM = 252G(y,0) - -F,{y) , 



Fi{y) = -2F2{v)., 



(Bla) 

(Bib) 

(Blc) 
(Bid) 



F^ = 



1 



F^{y) - -diG{y,H) - My) 



(Ble) 
(Blf) 



with s = 7ry/{2H). 
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2. Neumann boundary conditions 

For Neumann boundary conditions the functions appearing in the kernels in Eqs. ( pO[ ) and defined in Eqs. ( p2| ) have 
the exphcit forms 

T,{y) = 2G(y,0) + -1- £^ , (B2a) 



AnHy sinh(s) ' 

inh(s) 
16i?3y cosh3(s) ' 

TT CC 

16^^nh3(s) ' 
ih(s) 

AttHu cosh(s) 



T / \ ^ sinh(s) 



^ , N / N 7'' cosh(s) _ ^ 



^ . N 1 sinh(s) ,„ 
My) = rU, (B2d) 



^ / X 1 . ■"■'^ cosh^(s) TT^ cosh(s) , 

My) = --,ohiy,o) + + 161?^^^' 

, , TT^ sinh^(s) TT'^ sinh(s) , 

My) - - dloiy, ^) - ^ ^ + ^ ' 

^7(2;) = -G{y,0) + (B2g) 

Ts{y) = -G(2/,i/) + i^4(y), (B2h) 

My) = -2^7(y), (B2i) 

Miy) = -i^4(?/), (B2j) 

-^11(2/) = -^My), {B2k) 

^12(2/) - -9{y,H) - i^2(2/), (B21) 



with g{y,z) = diG{y,z) 
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